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We study the entanglement spectrum of highly excited eigenstates of two known models that 
exhibit a many-body localization transition, namely the one-dimensional random-field Heisenberg 
model and the qnantum random energy model. Our results indicate that the entanglement spectrum 
shows a “two-component” structure: a universal part that is associated with random matrix theory, 
and a nonuniversal part that is model dependent. The nonuniversal part manifests the deviation 
of the highly excited eigenstate from a true random state even in the thermalized phase where the 
eigenstate thermalization hypothesis holds. The fraction of the spectrum containing the universal 
part decreases as one approaches the critical point and vanishes in the localized phase in the ther¬ 
modynamic limit. We use the universal part fraction to construct an order parameter for measuring 
the degree of randomness of a generic highly excited state, which is also a promising candidate 
for studying the many-body localization transition. Two toy models based on Rokhsar-Kivelson 
type wave functions are constructed and their entanglement spectra are shown to exhibit the same 
structure. 
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Introduction.-Ulna.'DinTii entanglement, a topic of much 
importance in quantum information theory, has also 
gained relevance in quantum many-body physics in the 
past few years mm- In particular, the entanglement 
entropy provides a wealth of information about physical 
states, including novel ways to classify states of matter 
that do not have a local order parameter [3]. However, it 
has been realized only recently in various physical con¬ 
texts that the entanglement entropy is not enough to 
fully characterize a generic quantum state. For exam¬ 
ple, the quantum complexity corresponding to the geo¬ 
metric structure of black holes cannot be fully encoded 
just by the entanglement entropy [3]. One natural step 
beyond the amount of entanglement is the specific pat¬ 
tern of entanglement, i.e., the entanglement spectrum. 
A recent result that motivates this direction is the rela¬ 
tionship between irreversibility and entanglement spec¬ 
trum statistics in quantum circuits UM- It was shown 
that irreversible states display Wigner-Dyson statistics 
in the level spacing of entanglement eigenvalues, while 
reversible states show a deviation from Wigner-Dyson 
distributed entanglement levels and can be efficiently dis¬ 
entangled. 

Are there universal features in the entanglement spec¬ 
trum of a generic eigenstate of a quantum Hamiltonian? 
Highly excited eigenstates of a generic quantum Hamil¬ 
tonian are believed to satisfy the “eigenstate thermal¬ 
ization hypothesis” (ETH) [THH], which states that the 
expectation value {ifalOlipa) of a few-body observable O 
in an energy eigenstate \ipa) of the Hamiltonian with en¬ 
ergy Ea equals the microcanonical average at the mean 
energy Ea- So one could as well ask the following ques¬ 
tion: What is the structure of the entanglement spectrum 


of highly excited eigenstates of a thermalized system? 
Here we find a quandary. Completely random states 
are generically not physical, namely, they cannot be the 
eigenstates of Hamiltonians with local interactions. For 
the ETH to be a physical scenario for thermalization, 
highly excited eigenstates of physical local Hamiltonians 
cannot always be completely random, yet they have to 
contain enough entropy. Deviations from a completely 
random state can be quantified by the entanglement en¬ 
tropy, more precisely by the amount that it deviates from 
the maximal entropy in the subsystem, derived by Page, 
which we will refer to as the Page entropy hereafter m- 
But are there features that cannot be captured by the 
entanglement entropy alone? Can one identify remnants 
of randomness in the full entanglement spectrum? What 
about in states that violate the ETH? 

In this Letter, we address the above questions using 
as a case study the problem of many-body localization 
(MBL) [nHIS]- We study two known models that were 
shown to exhibit a MBL transition, namely, the Heisen¬ 
berg spin model with random fields, and the quantum 
random energy model (QREM) [T5II18) . In the delocal¬ 
ized phase, high-energy eigenstates are thermalized ac¬ 
cording to the ETH. The deviation from completely ran¬ 
dom states manifests itself in a “two-component” struc¬ 
ture in the entanglement spectrum: a universal part that 
corresponds to random matrix theory [TO], and a nonuni¬ 
versal part that is model dependent. We show that 
the universal part fraction decreases as one approaches 
the transition point and vanishes in the localized phase 
in the thermodynamic limit. We therefore propose an 
order parameter that is able to measure the degree of 
randomness of a generic highly excited state and cap- 
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ture the many-body localization-delocalization transition 
based on the entanglement spectrum, and show that it 
gives predictions consistent with previous results. We 
further construct two toy models in terms of Rokhsar- 
Kivelson- (RK) type wavefunctions [201111] and the same 
structure in the entanglement spectra is observed. 

Heisenberg spin chain.-A well-studied model that 
shows a MBL transition is the isotropic Heisenberg spin- 
1/2 chain with random helds along a fixed direction, 

L 

n = St + J 5, • 5,+i + ), (1) 

i=l 

where the random helds hi are independent random vari¬ 
ables at each site, drawn from a uniform distribution 
in the interval [—h,h]. F is a uniform transverse held 
along the x direction, which breaks total conserva¬ 
tion. We assume periodic boundary condition and set 
the coupling J = 1 and F = 0.1. In the absence of 
the transverse held F, previous work located the critical 
point dX h = he ~ 3.5 in the = 0 sector [I3l|22l|23]. 
We consider two different regimes by varying the disorder 
strength parameter h: (i) within the thermalized phase 
{h < he)-, and (ii) in the localized phase {h > he). In 
each regime, we focus on eigenstates of Hamiltonian Q 
at the middle of the spectrum, namely, on highly excited 
states. 

We consider a bipartition of the system into subsys¬ 
tems A and B of equal size {L/2 sites each). For a generic 
eigenstate |^) = where cr = (TiCT 2 ...o'l 

labels the 2^ possible spin conhgurations of the sys¬ 
tem, we cast the wave function as it(cr) = ipi^ACB), 
where cr^ = CTi... 0-^/2 and ctb = crB/ 2+1 ■ ■ ■ <^l- The 
entanglement spectrum is obtained from the eigenval¬ 
ues of the reduced density matrices pA = trB|'0)('0l and 
PB =tTA\i’){i’\-. {pfc = A|},A: = I,...2-^/2. In this work, 
we are primarily concerned with the density of states and 
level statistics of the {A^} for highly excited eigenstates 
for different strengths of disorder. For each value of h 
analyzed, the spectra were averaged over 10 realizations 
of disorder for L = 16, and 100 realizations for L = 14. 
For each spectrum, the eigenstate with energy closest to 
zero was obtained by a Lanezos projection [22]. This 
eigenstate corresponds to a highly excited state. 

Thermalized phase.-We start by considering the 
weakly disordered case, h he- Only a small amount 
of disorder is necessary to break the integrability of the 
clean Hamiltonian. However, conservation of the total 
also plays a crucial role in making eigenstates com¬ 
pletely random. A small transverse field F is applied 
to break this conservation without substantially altering 
the many-body localization transition. In this regime, 
we find that the entanglement spectrum of the highly ex¬ 
cited state with eigenenergy near zero is close to that of a 
completely random quantum state, as shown in Fig. Ba) 
for systems of size L = 16 and h = 0.5. The entangle¬ 


ment spectrum follows closely a Marchenko-Pastur distri¬ 
bution (with proper normalization), which describes the 
asymptotic average density of eigenvalues of a Wishart 
matrix [23 [2S|- (The expression for the entanglement 
spectral density for the random state is presented in the 
Supplemental Material |27|.l One can also check that, 
in this regime, the von Neumann entanglement entropy 
5(1)=-^ f.Pk Inpfc is in good agreement with the Page 
entropy for random states: Sm,n = YJk=n+i i ~ ~ 

ln(m) — where m and n are the Hilbert space di¬ 
mensions of subsystem A and B, respectively |10| . For 
example, our computed average entropy for 16 sites is 
= 4.9719 ± 0.0015, while the corresponding Page 
entropy is S'page = 5.0452. 

As the disorder strength is increased, but still h < he, 
the system remains in the thermalized phase where it is 
supposed to obey the ETH and yield volume-law scal¬ 
ing of the entanglement entropy with system sizes |30| . 
which is verified in the insets of Figs. Ba) to[^e). How¬ 
ever, in spite of the volume-law scaling of the entangle¬ 
ment entropy and the thermalization of eigenstates, the 
entanglement entropy is much lower than the Page en¬ 
tropy. This indicates that the pattern of entanglement 
must have changed, which is manifest in the spectra 
shown in Figs. Bb) to [^e). The entanglement spec¬ 
trum shows a striking “two-component” structure: (i) 
a universal tail in agreement with random matrix the¬ 
ory, and (ii) a nonuniversal part. The non-universal part 
dominates the weights in the spectrum (large Ak values), 
resulting in low entanglement entropy, as it decays much 
faster than the universal part. Therefore, we find that 
although thermalized states are not necessarily random 
states, they partially retain a component that is remi¬ 
niscent of a random state: the entanglement spectrum 
follows the Marchenko-Pastur level density distribution. 
In addition, the universal part of the entanglement spec¬ 
trum follows a Wigner-Dyson distribution of level spac- 
ings (see the Supplemental Material [2Z]). 

Localized phase.-hi this regime, the entanglement en¬ 
tropy exhibits an area-law scaling with the system size 
[see inset of Fig. , which in one spatial dimension im¬ 
plies a constant entropy and, at most, weakly logarithmic 
corrections, in accordance with Ref. [5T| . 

The entanglement spectrum in the localized regime, 
depicted in Fig. Bf) for h = 6, shows a different scenario 
from that in the thermalized phase: the universal part 
of the spectrum disappears completely, leaving only the 
nonuniversal part characterized by its fast decay rate. 

QREM.-The QREM describes L spins in a transverse 
field F with the following Hamiltonian: 

L 

H = E({a^})+F^a/ (2) 

i=l 

where E{{a^'\) is the classical REM term that takes inde¬ 
pendent values from a Gaussian distribution of zero mean 
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Figure 1: (Color online) Average entanglement spectrum of highly excited eigenstates for a system of size L — 16, averaged 
over 10 realizations of disorder (plotted in logarithmic scale). Panels a-f show the spectrum for h = 0.5,1.5, 2, 2.5, 3 and 6, 
respectively. The solid lines correspond to the spectrum of a completely random state (derived from a Marchenko-Pastur 
distribution), and is shown for reference. Insets: scaling of the average entanglement entropy with system size. 


and variance L/2 [32] . This model was first studied in the 
context of a mean-field spin glass, and was shown to ex¬ 
hibit a first-order quantum phase transition as a function 
of r [TB] . More recently, it was further demonstrated to 
have a MBL transition when viewed as a closed quantum 
system ffU- Numerical and analytical arguments show 
that the transition happens at an energy density |e| = T 
in the microcanonical ensemble. Since there is no sup¬ 
port for the many-body localized phase at energy density 
e = 0, we examine the eigenstates with energy density 
closest to e = 0.5 instead, and study the entanglement 
spectrum as T is tuned. The two-component structure 
and its evolution as a function of T similar to Fig. [Tjare 
again observed (see the Supplemental Material P7|). 

An order parameter.-The above picture unveils a new 
aspect of the MBL transition. The two parts of the entan¬ 
glement spectrum of a highly excited state clearly evolve 
as the disorder strength h is increased, namely, the uni¬ 
versal part shrinks and the nonuniversal part grows. This 
fact suggests that one could use the fraction of each com¬ 
ponent as an order parameter. 

Figures a) to[^ e) indicate an h dependent value kh 
that separates the nonuniversal (fc < kh) from the uni¬ 
versal (fc > kh) parts of the rank-ordered entanglement 
levels (see the Supplemental Material m for the proto¬ 
col for determining kh). One can thus define the partial 
Renyi entropies 

^ k<kh 

with 9 > 0. Because the universal part of the spectrum is 
where the eigenvalues with low entanglement reside, this 


part of the spectrum is obscured by any measure that 
relies on the eigenvalues as weights. A good measure of 
the fraction of the two components that does not depend 
on these weights is given by the <7 = 0 Renyi entropy, 
which simply measures the ranks: = Infc/j. There¬ 

fore, an order parameter that measures the fraction of 
the universal component is 


Ombl = 1 - = 1 - 


log2 kh 
L/2 


(4) 


Figure shows the order parameter as defined above 
for the Heisenberg spin model and the QREM, respec¬ 
tively. For the QREM, all curves at different system 
sizes cross at F^ « 0.5, in excellent agreement with Ref. 
m We have also looked at energy density e = 0.3, and 
the curves cross at Fc « 0.25, giving the same numer¬ 
ical prediction as in Refs. El and |18j (plot shown in 
the Supplemental Material [27|). For the random-field 
Heisenberg model, however, the fact that the transition 
happens at the point where the order parameter is nearly 
zero makes it harder to accurately locate the critical point 
using our order parameter. We see from Fig. [^that the 
curves cross at he ~ 3.3, which is also consistent with 
previous studies. This indicates that, by considering the 
full entanglement spectrum at high energies, our order 
parameter reveals a novel property that is promising for 
studying the MBL transition. 

We remark that, although the MBL transition can 
also be captured by the scaling property of the entangle¬ 
ment entropy, our order parameter seems to be applicable 
even for models with nonlocal interactions, which could 
obscure the connection between the volume-to-area law 
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Figure 2: (Color online) The order parameter defined as the 
fraction of the universal component in the full entanglement 
spectrum for the Heisenberg spin model (upper panel) and 
the QREM (lower panel). 

transition of the entropy and the MBL transition. 

Toy models.-we construct two RK-type model wave 
functions that are shown to have (i) the two-component 
structure in their entanglement spectra, and (ii) a phase 
transition as a function of the tuning parameter. The 
wave functions take the following form: 

= ( 5 ) 

where E{(t) is the energy for the classical configuration cr 
and Z is the corresponding partition function of the clas¬ 
sical statistical system [21] . 5^. is a random sign for each 
configuration, such that the wave function represents a 
highly excited state. We consider the following two cases: 
(i) E{(t) = E^em{(t), and (ii) E{(t) = -{ T,i<j ■ 
In the first case, the energy is taken to be that of the 
REM, while in the second case the energy is that of an 
infinite-range uniform ferromagnetic interaction. 

In the small /3 regime, the above RK-type wave func¬ 
tions are close to completely random states; upon in¬ 
creasing /3, the wave functions are pushed towards prod¬ 
uct states and start to deviate from completely random 
states. Therefore, the tuning parameter /3 here plays the 
role of the “disorder strength”. Indeed, we find the same 
two-component structure in the entanglement spectrum 
(see the Supplemental Material [17]), and the order pa¬ 
rameter is shown in Fig. [^ The REM case was recently 


studied by Chen et al. where the MBL transition was 
obtained numerically using other measures [33] . Here we 
clearly see that, in both cases, the curves cross at some 
critical [3, indicating the existence of a similar phase tran¬ 
sition. 

Summary and discussion.-The details of the structure 
of the entanglement spectrum, especially the universal 
part at the tails of the spectrum, have long been over¬ 
looked. The main focus has been primarily on the dom¬ 
inating nonuniversal component, and the universal tail 
has thus far been discarded. For example, in the den¬ 
sity matrix renormalization group [34] and tensor net¬ 
work methods [35] . the density matrix is truncated to 
avoid uncontrolled growth of its dimensions. While this 
procedure is certainly justified when the purpose is to 
obtain ground state properties, it discards important in¬ 
formation about the behavior of the system at higher 
energy states. In this Letter we showed that the full en¬ 
tanglement spectrum, directly computable from the wave 
function, provides information that is often invisible in 
the entanglement entropy alone. 

On the other hand, much has been known about ran¬ 
dom quantum states, e.g., the Page entropy and volume- 
law scaling entropy. Nevertheless, the Page entropy is 
often an overestimate of the actual entanglement entropy 
computed from generic quantum states. Therefore, a nat¬ 
ural question that arises is as follows: How random does 
a given quantum state look? In this Letter, we show 
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Figure 3: (Color online) The order parameter for the random- 
sign RK-type wave functions. Upper panel: E((t) = 
^Irbm(o'). Lower panel: E{(t) = - ^ J = 1. 
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that a generic quantum state that satisfies ETH does not 
necessarily mean a completely random state. We present 
an order parameter to quantify the degree of randomness 
by using information about the full entanglement spec¬ 
trum. In the context of MBL, our order parameter is 
able to locate the critical point, consistent with previous 
results. Our work may provide a novel way of study¬ 
ing MBL, and may shed new light on the understanding 
many-body systems at the level of wave functions. 
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SUPPLEMENTAL MATERIAL 

MARCHENKO-PASTUR DISTRIBUTION AND RANDOM STATES 


The Marchenko-Pastur distribution describes the asymptotic (large-A) average density of eigenvalues of an A x A 
matrix of the form Y = AA^, known as a Wishart matrix, where A is a A x M random rectangular matrix with 
independent but identically distributed entries [25]. Let cr^ be the variance of the entries in A. When N = M ^ oo, 
the Marchenko-Pastur distribution takes the form 


\ k=l / N^ao * 

where {pk} are the eigenvalues of A, 0 < p < Pmax and Pmax = 2cr^. From this distribution we can obtain the average 
number function associated to the eigenvalues of Y. Let pi > P 2 > • • • ^ Pn and rjk = k/N. Then, 


vip) = 1 

= 1- / dp'V{p') 

^0 


N—¥00 


dp'V{p') 


Jp 

TT L 


u\/1 — -I- arcsin(u) 

Thus, introducting the rescaled variable x = y/p/pmux, we hnd 


“ = -\/p/Pn 


r]{x) = I- x\/l — -I- arcsin(a:) 


(7) 

( 8 ) 
(9) 

( 10 ) 


( 11 ) 


It is straightforward to relate the average number function derived from the Marchenko-Pastur distribution with 
that obtained from the entanglement spectrum of a bipartitioned random vector. Let iI){xa,Xb) be the wavefunction 
of the bipartite system. Then, the reduced density matrix is given by 


Pa{xa,x'a) = y^^1p{xA,XB)'tp*{x'A,XB)- 

Xb 


( 12 ) 


We can see that, for completely random wavefunctions, the reduced density matrix is a random Wishart matrix and 
therefore its eigenvalues should follow a Marchenko-Pastur distribution[26]. Thus, we expect the average number 
function to provide an accurate description of the average spectrum. 

Let {pk}, k = l,...,d, be the set of eigenvalues of pA in decreasing order, with pk > 0, J2k=iPk ~ 1> ^'^‘1 
d < 2^/^. It is straightforward to relate the eigenvalues {pk} to the singular values {Afe} resulting from the Schmidt 
decomposition of the bipartite wavefunction. 


-IpixA.XB) = 


d 

E 


Afc (j)A\^A)(pP{xB), 


(13) 


by simply setting pk = (notice that Xk >0), where 4>^^\xa) and (p^g’(xB) are the left-singular and right-singular 
vectors, respectively. For the purpose of comparing the average spectra to Eq. (II), it is necessary to rescale the 
singular values and their indices as follows: 


k(k), 



( 14 ) 
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where x{r]) is the inverse function of r]{x). The prefactor is chosen to guarantee the normalization of pk- 


d pX 

'^Pk ^ 4: dpx'^iT]) 

fc=i -^0 


= — / dxx^\/l —. 


TT 

= 1 . 


(15) 


We tested this formulatio n b y plotting the numerical results for pk obtained from a random state against the 
analytical expression in Eq. (111. Figure 4 shows ?7fc = g versus Xk = ^y^pkd. There is very good agreement with the 
analytical prediction. 



Figure 4: Average entanglement pectrum of completely random wavefunction (L = 16, 100 realizations used). 

Notice that in the cases considered in the main text, the random part of the entanglement spectrum alone is not 
normalized. However, by plotting logA^ versus fc, the missing normalization factor only amounts to a trivial shift of 
the entire spectrum. 


LEVEL SPACING STATISTICS 

We studied the statistics of the entanglement spectrum by looking at the level spacing distribution in the set {A*,}. 
To avoid having to perform spectral unfolding, we chose to evaluate the distribution of ratios of adjacent level spacings 
m- Xk = (Afc +1 — Afc)/(Afc — Afc_i). Accurate surmises exist for the distribution of these ratios in the case of Gaussian 
ensembles [2S]. They are given by 


, . 1 (r + 


(16) 


where Z = ^ iov the Gaussian Orthogonal Ensemble (GOE) with ,0 = 1, and Z = for the Gaussian Unitary 

Ensemble (GUE) with (3 = 2. The corresponding distribution for the Poisson distributed spectrum is given by the 









exact form 


1 


-f^Poisson — 


(1 + r )2 ■ 
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(17) 


Notice that for the Gaussian ensembles, level repulsion manifests itself in the asymptotic behavior P{r —)■ 0) ~ 
which is absent in the case of Poisson statistics. 

Results are shown in Fig. for disordered Heisenberg chains with L = 16 and 100 disorder realizations. A 
completely random real state follows a GOE statistics. The universal part of the spectrum at h = 0.5 also follows a 
GOE distribution. 




Figure 5: Left panel: Distribution of the ratios of consecutive spacings for the entanglement spectrum of Heisenberg spin chains 
with disorder parameter h = 0.5 and transverse field F = 0.1 (crosses). Right panel: the same data in a logarithmic scale. 
L = 16 and 500 realizations used. 


PROTOCOL FOR DETERMINING kn 


The definition of the order parameter in this Letter required a protocol for determining the point kh which separates 
the non-universal part (k < kh) from the universal (fc > kh) part of the rank-ordered entanglement levels. We use the 
following protocol: we took the spectrum obtained from each random state considered and multiply it by a factor s 
such that the rescaled smallest singular value coincided in value with that obtained from a completely random state. 
Then we swept through the spectrum, starting from the tail, and computed the relative deviation from the completely 
random state prediction, until it exceeded a certain amount. That is, until 




> e, 


(18) 


where e is a number of order 1. In our case, we set e = 1. 

However, we would like to point out one subtlety of this methodology. In cases where the universal component of 
the spectrum almost vanishes near the transition, it is hard to accurately locate the critical point. That is because the 
last few points at the tail of the spectrum show large sample-to-sample fluctuations, and our protocol requires strictly 
matching a single point close to the tail. Therefore, kh can be very sensitive to our choice of the matching point and 
can even yield incorrect predictions under finite-size scaling. On the other hand, we find that in cases such as the 
QREM, where there is still a large fraction of the universal component at the transition, this methodology is not very 
sensitive to the choice of matching point. For the Heisenberg model, we locate the critical point kh by choosing the 
matching point away from the tail end, thus effectively discarding the smallest eigenvalues. For example, for L = 16 
we discarded the last 16 eigenvalues. In order to demonstrate that this does not lead to a sizable errors in determining 
the MBL critical point, we also computed the the order parameter directly from the averaged entanglement spectrum. 
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where the fluctuations are smoothened out (Fig. [^. The critical point found this way is very close to the one we show 
in the main text. We believe that the result can be further improved if more realizations are include in the averaring, 
which we hope to attempt in the future. 





Figure 6: Left: Results obtained when the matching point is chosen too close to the last eigenvalue for each entanglement 
spectrum, yielding a wrong transition point under finite-size scaling. Middle: the order parameter computed from the averaged 
entanglement spectrum of the random-field Heisenberg model, where fluctuations are smoothened out and not many small 
eigenvalues need to be eliminated. Right: the order parameter computed from the average entanglement spectrum of the 
QREM. One can see that in this case our approach is robust and the result is very close to that shown in the main text. 


ENTANGLEMENT SPECTRA FOR QREM AND RK-TYPE TOY MODELS 

In this section, we present in Fig. [^through Fig. [^the entanglement spectra for the QREM and two RK-type toy 
models that were discussed in the main text. 



Figure 7: Averaged entanglement spectrum of eigenstates with energy density e = 0.5, for the QREM. The system size L = 16, 
and r = 0.9, 0.6, 0.4, and 0.2, averaging over 10 realizations of disorder. 

We clearly see the following: (1) the two-component structure shows up in all three cases; (2) the same evolution 
behavior as explained in the main text happens here as well. Namely, the universal fraction shrinks as one increases 
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Figure 8: Averaged entanglement spectrum of the RK-type wavefunction, with E{cr) = -Erem(o'). The system size L = 16, 
and P = 0.2, 0.7, 1, and 1.8, averaging over 100 realizations of disorder. 





Figure 9: Averaged entanglement spectrum of the RK-type wavefunction, with E{cr) = — ^ The system size 

L = 16, and /? = 0.3, 1, 1.5, and 2.5, averaging over 100 realizations of random sign. 


the strength of disorder, thereby pushing the states further away from completely random states. 

We also show the order parameter for the QREM at energy density e = 0.3, which is different from the one shown 
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in the main text. We clearly see from Fig. 10 that the curves for different system sizes cross at around Fc « 0.25, 
which is again in excellent agreement with previous known results. 



Figure 10: The order parameter for QREM with target energy density e = 0.3. 


RANDOMNESS VERSUS NON-RANDOMNESS: ANOTHER TOY MODEL 

In this section, we view the emergence of the non-universal component in the entanglement spectrum from a different 
perspective: the degree of randomness in the wavefunctions. The physical intuition can be understood as follows. 
A completely random state is supposed to yield a Marchenko-Pastur distribution in its entanglement spectrum, i.e. 
only the universal component exits. This implies that for states whose entanglement spectra deviate from Marchenko- 
Pastur distribution, they cannot be completely random. Therefore we construct another toy model which captures 
this feature, by borrowing ideas from spin glasses. 

First let us start with a truly random (real) wavefunction that we denote by = sgn(AREM)j where sgn(a;) 

is the sign function and Arem are identically independently distributed with probability P(£'rem) = 
with L being the number of spins. The wavefunction T'rem iS; by construction, random. The subscript REM is used 
to draw an analogy to the Random Energy Model (REM) in spin glass systems [? ]. 

Next, we note that the REM is a limiting case of a spin glass with p-spin interactions as p —> oo, which elim¬ 
inates correlations between configurations. So we can take a step back and consider the following “less random” 
Sherrington-Kirkpatrick (SK) spin glass model with infinite-range two-spin interactions |29j . and construct the fol¬ 
lowing wavefunction: 


^SK(cr) = sgn(E;sK(cr)) = sgn ^ J,. 


(19) 


, i<3 


where the are drawn from uniform distribution in the interval [—1,1]. The amplitudes computed from the SK-like 
model are obviously not as random as in the REM-like one; there are only L(L — l)/2 independent random J^’s 
in the former as opposed to 2^ independent random amplitudes in the latter. Nevertheless, the amplitudes of T'sk 
do inherit some randomness from the and the energy distribution of the SK-like model also follows accurately a 
Gaussian distribution for any one given cr, similarly to those in the REM-like model. But the correlations between 
the amplitudes for different cr exist in the case of the SK-like model, and these correlations are manifest in the 
entanglement spectrum computed from T'sk, as shown in the left panel of Fig. [3 

The entanglement entropy follows a volume-law scaling (see inset of Fig. 11) and, again, we see the emergence of a 
two-component structure in the spectrum. The universal part agrees with RMT. The non-universal part is different 
from that found for the high energy eigenstates of a disordered Heisenberg spin chain, reflecting its non-universal, 
model dependent nature. Yet, this component is still characterized by its fast decay rate. The toy model shows that 
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Figure 11: (Color online) Entanglement spectrum of: ^sk (left) and ^ 3 _apin(fT) (right), for a system of size L = 16, averaged 
over 500 realizations of disorder. The order parameter as defined in Eq. (3) is: {O) = 0.3888 ± 0.0067, with a threshold k ~ 30 
for ^sk; and (O) = 0.4279 ± 0.0060, with a threshold k ~ 24 for 4'3_apin(o’). The insets show the volume-law scaling of the 
von Neumann entanglement entropy. 


non randomness can be present in a generic quantum state when there are correlations between components of the 
wavefunction. The entanglement spectrum captures the non randomness and its structure can be well described by a 
two-component picture. 

Another interesting manifestation of the mixing between universal and non-universal components in the SK wave- 
function is revealed by employing a color map. In Fig. 12 we show the amplitude of the wavefunction 1 E'sk(c’’a, (^b) 
plotted in a cr^ x ctb grid, and compared it to the amplitude of a REM wavefunction. The existence of a struc¬ 
ture, similar to wefts in a tapestry, is clearly visible for the T'sk wavefunction, but completely absent for the REM 
wavefunction. 



Figure 12: (Color online) Color map of the matrix ’1 '((Ta, (Tb) for typical realizations of the (a) SK and (b) REM wavefunctions 
(L = 18). 

One can further consider wavefunctions built with three-spin interactions, 


4'3_spi„(cr) = sgn ^ Jijk CTi (7j Ok | , 

\i<j<k 

where the Jijk are drawn from a uniform distribution in the interval [—1,1]. As shown in the right panel of Fig. 
the spectrum again shows a two-component structure, very similar to the cases discussed in the main text. 
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